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ABSTRACT 

The two dominant features in the distribution of orbital parameters for close-in exo- 
planets are the prevalence of circular orbits for very short periods, and the observation 
that planets on closer orbits tend to be heavier. The first feature is interpreted as a 
signature of tidal evolution, while the origin of the second, a "mass-period relation" for 
hot Jupiters, is not understood. In this paper we re-consider the ensemble properties 
of transiting exoplanets with well-measured parameters, focussing on orbital eccen- 
tricity and the mass-period relation. We recalculate the constraints on eccentricity in 
a homogeneous way, using new radial-velocity data, with particular attention to sta- 
tistical biases. We find that planets on circular orbits gather in a well-defined region 
of the mass-period plane, close to the minimum period for any given mass. Excep- 
tions to this pattern reported in the Literature can be attributed to statistical biases. 
The ensemble data is compatible with classical tide theory with orbital circularisation 
caused by tides raised on the planet, and suggest that tidal circularisation and the 
stopping mechanisms for close-in planets are closely related to each other. The posi- 
tion mass-period relation is compatible with a relation between a planet's Hill radius 
and its present orbit. 
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1 INTRODUCTION 

Transiting planets are presently our main source of infor- 
mation on the formation, structure and evolution of extra- 
solar planet^] Thanks to the success of photometric transit 
surveys, enough transiting planets are now known to start 
applying more robust statistics to their ensemble features. 
In particular, transiting planets help to address the question 
of two puzzling features of the presently known population 
of exoplanets: the presence of planets on very close orbits, 
and the wide range of orbital eccentricities. 

Because of the strong selection biases favouring shorter 
orbital periods, most known transiting planets have close 
orbits (a < 0.1 AU). Contrary to the sample of planets 
found by radial-velocity surveys at larger orbital distances, 
the orbits of transiting planets tend to be more circular, a 
tendency that is interpreted as a signature of tidal circu- 



larization (Mazeh 20081. A close study of the eccentricity 



1 see compilations in www.exoplanets.eu (unedited catalogue of 
published parameters), www.exoplanets.org (edited catalogue fo- 
cussing on Doppler detections), www.inscience.ch/transits (edited 
catalogue of transiting planets), and in Southworth I 2009 1. 



distribution of these objects needs to address two related 
issues: how did close-in planets get there, and what is the 
role of tidal orbital evolution? Early inward migration by 
gravitational interaction with the protoplanetary gas disc is 
the favoured explanation for the presence of close-in planets 
( |Lin et al.|1996[ ), with planet-planet scattering, followed by 
tidal dissipation, as another possibility | Rasio fc Ford|1996 1 . 

The empirical evidence that has to be accounted for 
by a successful understanding includes the pile-up of close- 
in planets near orbital periods of 3 days, the predominance 
of circular orbits inwards of a few days, and the possible 
existence of a correlation between mass and period for close- 
in gas giants (with more massive hot Jupiters found closer 
in, P—l-2 days rather than 3-4 days), as pointed out by 



Mazeh et al. (2005) 



The classical theory of tides predicts that the orbits of 
transiting planets will be circularised by the tidal interaction 
induced by the star on the planet, with a timescale increas- 
ing sharply with orbital distance (e.g. |Goldreich fc Soter| 
|1966[ ) . The observed transition from eccentric to circular or- 
bits around periods of a few days is qualitatively compatible 
with this prediction. Similar behaviour has been found for 



binaries, e.g. Mathieu & Mazeh (1988) and Mazeh (2008) 
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Figure 1. Period-eccentricity plot for transiting planets using a 
compilation of the latest values published in the literature (top) 
compared to the determinations in the present study (bottom). 
The 95% lower intervals are shown when a null eccentricity is not 
excluded at the 5% level. Note the difference for low-eccentricity, 
short-period values. 



Recently though, a number of small but non-zero or- 
bital eccentricities have been reported for several planets on 
extremely close orbits (below 3 days), including WASP-12b 
( Hebb et al.|[2009b, WA SP-14b ( |Joshi et al |[2009|), WA SP- 
18b ( [Hellier et al.|2009[), WA SP-19b ( |Hebb et al.|2010[ ) and 
GJ 436b ( |Butler et al.|2004[ ) — see top panel of Fig. |lj Spe- 
cific scenarios have been invoked to explain these apparently 
anomalous cases, such as the presence of a perturbing com- 
panion ( |Ribas et aL]|2008| jBatygin et al.||2009[ ), or widely 
different coefficients in the intrinsic response of planets to 
tides (Matsumura et al. 112009 1. Another relevant recent de- 



velopment is the detection of several systems with strong 
spin-orbit misalignment (Winn et al. 2010 and references 



therein), which is presently interpreted as an indication that 
the disc-migration scenario is in need of significant updates, 
and that dynamical evolution after the dissipation of the 
disc is probably more important than previously thought, 
and even possibly dominant in determining the orbital prop- 
erties of close-in planets. 

We have been performing a series of observational pro- 
grammes to gather more radial-velocity (RV) data on known 
transiting planets, in order to study these issues. In this pa- 
per, we revisit the data on orbital eccentricity for transiting 
planets, based on new data and a re-analysis of extant data, 
and examine the implications in terms of the ensemble prop- 
erties of close-in planets regarding the orbital circularisation 
and the stopping mechanism. 



2 PLANET SAMPLE AND ORBIT ANALYSIS 

Our sample consists of all known transiting planets fulfill- 
ing the following selection criteria: (i) objects with precisely 
determined parameters (i.e. planetary radius and planetary 
mass measured to better than 10% accuracy) (ii) host stars 
brighter than V=15 mag (iii) objects known from peer- 
reviewed publication by 1st July 2010. 

We use published radial velocity data as well as new 



HARPS data from our observation programme for WASP- 
2, WASP-4, WASP-5 and WASP-7. The radial velocity data 
for these objects is presented and discussed in more details 
in Husnoo et al. (20111, and is available as an electronic 



table. 

The brightest known transiting planets were identified 
by Doppler planet surveys first and only subsequently found 
to be transiting. Their orbital eccentricity is therefore usu- 
ally well determined, because a sizeable number of individual 
radial velocity measurements is required to find the period 
of a planetary orbit. Most of our sample, however, consists 
of objects discovered by photometric searches for transits. 
In that case, the number of radial velocity measurement is 
lower, since the RV data are only required to detect the 
presence of an orbital signal at the period and phase of 
the transits. Transiting planets from photometric searches 
are on average much fainter than the targets of Doppler 
searches, so that in general the RV data is also less precise. 
In these case, there are three options to calculate the orbital 
parameters: 

- the eccentricity is set to zero, under the assumption that 
it is small and undetectable, and the data is solved for the 
other parameters of the system (e.g. the planets from the 
OGLE, HAT and TrES surveys) 

- the RV data clearly indicate an eccentric orbit, and the 
eccentricity e is solved together with the other parameters 
(e.g. HAT-P-2). 

- a Monte-Carlo Markov Chain (MCMC) integration is run 
on all data, with the eccentricity as a free parameters (e.g. 
planets from the WASP survey). 

The problem with the third method is that there is 
an intrinsic bias in the determination of eccentricity from 
radial-velocity data, because measurement errors on RV 
data for a circular orbit will always give rise to a finite 
best-fit eccentricity. This effect was studied in the context 



of stellar binaries by Lucy & Sweeney ( 1971 1 and more re- 



cently by Shen & Turner (20081. As a result, even when 



the measurement uncertainties are correctly estimated, the 
centre of the posterior distribution from a MCMC for a cir- 
cular orbit will tend to be 1-3 sigma away from zero. An 
additional issue is the presence of other sources of signal in 
the radial velocity data. The most common causes are stel- 
lar activity, instrumental drifts, and additional companions 
in the system. When these sources are not included in the 
orbital analysis, they may induce spurious eccentricity de- 
tections (Rodigas & Hinz 2009" Husnoo et al.|20lT | because 
any radial-velocity offset from a circular orbit will make it 
appear more eccentric. 

A large fraction of published eccentricities for transiting 
planets are significantly different from zero at the few-sigma 
level only. Particular attention to this issue is justified to 
avoid spurious detections, which would be especially con- 
fusing in the context of studying tidal evolution, since they 
would tag some circularised objects as having low but sig- 
nificant eccentricity, thus being good candidates for rapid 
on-going tidal evolution. 

The recent case of WASP-12, di scussed in|Hebb et al.| 
([20091); |L6pez-Morales eTaL] poTof ; |Campo et al.| ( |2010[ ); 



Husnoo et al. (20111, shows how false positives can arise 



in non-circular orbit detections. The orbital eccentricity of 
WASP-12b, initially thought to be significant at the 3-sigma 
level, was subsequently found to be spurious by further RV 
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measurements and the detection of the secondary eclipse by 
the Spitzer satellite. Another example is WASP-10, recently 



re-analysed by Maciejewski et al. (20101. These authors con- 
clude that the initially significant eccentricity detection was 
in reality due to radial- velocity variations induced by stellar 
activity. 

For these reasons, we re-analyse the RV data for tran- 
siting planets with a specific attention to the discrimination 
between circular and eccentric orbits. 

We use a Markov Chain Monte Carlo (MCMC) method 
with the Metropolis-Hastings algorithm (see Collier- 



Cameron et al. 2007 for example). Our implementation is 
described in Pont et al. (20091. We take into account the 



possible non-random RV uncertainties in addition to the for- 
mal error bars, from either physical or instrumental causes, 
ensuring that the size of the O — C residuals is compati- 
ble with the uncertainties used for the MCMC. We include 
the possibility of a certain amount of correlation in the ad- 
ditional RV noise (using the methods in Pont et al.| [2006). 
Finally, to decide whether the evidence for an eccentric orbit 
is statistically significant, we examine the "e = 0" possibil- 
ity in hypothesis-testing mode rather than in parameter- 
estimation mode. In other words, we do not simply adopt 
the best-fitting value of e as the most likely value, but also 
calculate if a non-zero value of e is more likely than the 
e — null hypothesis given the addition of free parame- 
ters and the possibility of non-random noise in the data. We 
do this by calculating the BIC (Bayesian Information Crite- 
rion) of the best-fit eccentric orbit compared to the best-fit 
circular orbit. An eccentric Keplerian has two extra free pa- 
rameters compared to a circular one, so that a closer fit is 
always expected even in the absence of real evidence. This 
last step penalises the eccentric solution for the two extra 
parameters. 

In summary, we consider the evidence for orbital eccen- 
tricity to be sufficient not only when an eccentric orbit gives 
a better fit to the data at a level higher than the uncer- 
tainties, but when the BIC value shows that the improve- 
ment compared to a circular orbit is significantly larger than 
would be expected from instrumental and physical sources 
of noise and added free parameters. 

To circumvent the difficulty posed by the tight correla- 
tion of the e and ui parameters near e = 0, we use ecosw 
and esincj as parameters instead of e and uj, following |Ford| 
(20061. Using these parameters also makes it easier to avoid 



biasing results towards higher eccentricities when expressing 
the output of the MCMC chain in terms of confidence inter- 
vals. The posterior probability distribution for e is one-sided 
near because e < is forbidden, but e cos uj and e sin uj do 
not have this problem. We calculate the central values of the 
probability distribution of e and uj from the central values 
for ecoscj and esinu; in the output of the MCMC. 



3 RESULTS: ORBITAL ECCENTRICITIES 
FOR TRANSITING PLANETS 

Table [l] shows the new eccentricity determinations. When 
the detection of eccentricity is significant, the value deter- 
mined with the MCMC fit for the eccentric orbit is given, 
together with the 68% central confidence interval. For sys- 
tems where our analysis showed there was no evidence for 



an eccentric orbit, the upper 95% confidence interval on the 
eccentricity is given. 

Fig. [I] compares our eccentricity determinations with a 
compilation of published values. Some objects were already 
analysed with methods similar to our own, and in that case 
we adopt the eccentricity distribution from the Literature 
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published values, our new RV data and re-analysis result in 
better constraints of the eccentricity for several objects, and 
significant changes for others. But the main result is that 
in several cases, an eccentricity detection at the few-sigma 
level is transformed into an upper limit, i.e. the evidence for 
orbital eccentricity is found to be insignificant. Specifically, 
there are 10 transiting hot Jupiters with eccentricity detec- 
tions significant at more than the 2-sigma level below period 
of 5 days in the literature. In our re-analysis, there are five 
below 5 days, none of which are in the 0.5-2 Mj mass range 
typical of hot Jupiters. The systems with >2<r published ec- 
centricity detections that we do not find to be significant are 
CoRoT-5, HAT-P-13, WASP-5, WASP-6, WASP-10, WASP- 
12, WASP-17, WASP-18. 

There are two main causes for this: either new RV mea- 
surements show the initial eccentricity detection to have 
been spurious, or the hypothesis-testing statistic shows that 
the difference between the best-fitting eccentric orbit and 
a circular solution is not significantly larger than would be 
expected from intrinsic bias and measurements errors, and 
therefore a circular orbit is compatible with the data. With 
our new RV data, we have found evidence of RV modu- 
lations caused by stellar activity, instrumental effects or an 
additional companion in the system for WASP-4 and WASP- 
5. For objects measured with the SOPHIE spectrograph in 
the "HE" mode (used for fainter objects), we also find large 
long-term instrumental drifts, that need to be included in 
the error budget. 

The most striking feature of Table [I] is the scarcity of 
confirmed eccentricity detections below the 10% level com- 
pared to published values (see Fig. [I]). WASP-14, with e = 
0.088 ±0.003 and P = 2.2 days, becomes the shortest-period 
orbit with significant eccentricity detection, and HAT-P-16 
(e = 0.035 ± 0.003) the smallest detected eccentricity. 

Note that the scarcity of confirmed e < 0.1 orbits is not 
simply a reflection of the fact that smaller eccentricities are 
more difficult to detect. For several objects, the uncertainty 
on e is smaller than 1 percent, which would have allowed the 
detection of a small but non-zero eccentricity (this is the case 
either for bright and well-measured objects like HD 189733, 
or for objects with Spitzer secondary eclipse measurements 
like WASP-18). 

For lower-mass planets (below 0.5Mj) the eccentric- 
ity measurement gets difficult because the RV amplitude 
is smaller, and the constraints on e are rarely good enough 
to be significant in our context. That is a bias to keep in 
mind when studying the sample as a whole. We therefore 
need to compare significant eccentricities with eccentricities 
compatible with zero only when the uncertainties are small 
enough that a significant eccentricity would have been de- 
tected. 
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Figure 2. Mass-period diagram for transiting planets. Open circles (o) indicate objects with well-constrained orbital eccentricity com- 
patible with zero, filled circles (•) show definitely eccentric orbits with e > 0.1, and crosses (+) objects for which no definite conclusion 
could be derived from the present data. Gray circles ( ) indicate objects with a detected non-zero but small eccentricity (compatible 
with e < 0.1, values that could remain undetected at lower masses. The dotted line is the relation from|Mazeh et aJT] (|2005[l . 



4 DISCUSSION: TIDAL CIRCULARISATION 
AND THE MASS-PERIOD RELATION 

Figure [2] displays our sample in the mass-period diagram. 
The different symbols indicate probably circular orbits, ec- 
centric orbits, and unsolved cases. In the "probably circular" 
category, we include all objects with RV data compatible 
with a circular orbits and excluding an eccentricity larger 
than 0.10 at more than the 3-sigma level. 

Figure [2] shows a clear relation between orbital eccen- 
tricity and the position in the mass-period diagram, much 
clearer than is the case when a compilation of published 
values of eccentricities is used. In fact, there is no confirmed 
eccentric orbit at all in the main "hot Jupiter clump" (at 
P=l-4 days and M=0.5-2 Mj). A contrario, all objects that 
are on wider orbits and/or are more massive have eccentric 
orbits. 

Another conclusion is that mass and period are corre- 
lated for the objects with circular orbits — if not by a tight 
relation as proposed by |Mazeh et al.| ( |2005[ ), then definitely 
a trend, with heavier planets closer in. Especially remark- 



able is the tight grouping of lighter (M < Mj) planets on 
circular orbits. 

These two features are our main results: close-in planets 
belong to two clear classes in regard to orbital parameters - 
either they follow circular orbits and a mass-period trend, or 
they have eccentric orbit and are off the mass-period trend 
(on the side of larger period at a given mass). This is very 
strongly suggestive of tidal effects. 

Theoretical treatments of tidal evolution predict that 
the timescale of tidal circularisation will decrease with the 
planet-to-star mass ratio, and increase steeply with the semi- 
major axis scaled to the planetary radius. For the orbital cir- 
cularisation of planetary orbit, the relevant effect is the tide 
raised on the planet by the star rather than the opposite, see 
e.g. Hansen et al. 2 010 for a quantitative discussion. |Goldre^ 



ich & Soter 



|1966| obtain r e = ±Q (^L) £ 

tor the circularisation timescale, where Q is the tidal qual- 
ity factor — a parametrisation of the response of the 
planet's interior to tidal perturbation, r e is the circulari- 
sation timescale, G the gravitational constant and a the or- 
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Figure 3. "Tidal" diagram for transiting planets, showing the mass ratio of the system M p /M* versus orbital distance in terms of 
planetary radius (a/r p ). Same symbols as on Fig. [5] The dashed line shows a 1 Gyr circularization timescale isochrone, using Q = 10 6 , 
P = 3 days and e = 0. The dotted line shows a = 2 • an, the location predicted by tidal circularisation from very eccentric orbits (see 
text). 



bital semi-major axis. When the objects in our sample are 
ranked by T e /Q, we find that circular and eccentric orbits 
are separated as expected from classical tide theory, with an 
overlap in r e /Q between the smallest measured eccentrici- 
ties and the circular orbits over less than one decade, that 
can be accounted for by expected variations in Q and age 
from object to object and do not require additional factors 
like perturbing companions or anomalous tidal evolution. 

Figure [3] shows the position of our sample in a plane 
more closely related to the circularisation timescale than 
the mass-period diagram. Figure [3] displays the planet-to- 
star mass ratio M p /M*, as a function of the orbital distance 
scaled to the planetary radius a/R p . If R p and M a are much 
less variable than M v and P, and P and a are closely re- 
lated, as is the case for hot Jupiters, this plot will be similar 
to the mass-period plot, except that the tidal circularisa- 
tion isochrones will be straight lines. The dashed line shows 
a 1 Gyr circularisation isochrone using Q=10 6 (a plausible 
value for Jovian planets, Wu ( 2005 1), P — 3 days and e ~ 



(the dependence of the circularisation timescale on P and e 
is weak on such a log-log plolQ. On this plot, the relation 
between circularization timescale and measured eccentricity 
is apparent. Circular and eccentric orbits are well separated 
and the separation is compatible with circularisation on a 
timescale corresponding to the age of the systems. The small 
but non-zero eccentricities are all in the intermediate range. 
Thus, the whole sample is compatible with the operation of 



2 The timescale depends on e because an eccentric orbit has a 
lower angular momentum than a circular one at a given period. 



tidal circularisation with less than one decade of variation in 
system ages and values of Q, without any exception. This is 
very remarkable given the variety of planets and host stars 
considered. The slope of the relation shows that orbital cir- 
cularisation is caused by tides raised on the planet, rather 
than tides raised on the star that would predict a mass de- 
pendence in the opposite direction. 

The locus of the separation between circular and eccen- 
tric orbits is not compatible with timescales comparable to 
the lifetimes of protoplanetary discs (^ lOMyr). This sug- 
gests that the stopping mechanism for inward migration is 
not related to disc migration, or a coincidence would be re- 
quired to explain the correspondance of circular orbits and 
the mass-period relation. 

Another noteworthy element is the scarcity of high-mass 
planets on circular orbits ( Rasio et al.||1996 Jackson et al 



2009| | Winn et al. 20101. Although a larger sample would 
be necessary for definite conclusions, it seems that the re- 
sponse to orbital circularisation has a strong dependence on 
the mass of the planet: M<lMj planets stop at the circular- 
isation radius, M~lMj planets can move closer to the host 
star, and M>lMj planets do not reach circularisation. A 
possible explanation for this behaviour is that heavier plan- 
ets could raise tides in the star strong enough for angular 
momentum exchange and orbital decay ( Rasio fe Ford|1996 



| Sasselov|[2003l |Jackson et "alj]2009l |Pont||2009| |Winn et al.| 
|2010[ > . Another possibility is that heavier and lighter gas gi- 



ants form by different mechanisms (Ribas & Miralda-Escude 



2007 


I- 




Ford & 


i Rasio| 



Ford &: Rasio| ( |2006| > pointed out that if the orbit of 



close-in planets evolve from initially very eccentric orbits 
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due to interactions with other planets, the final orbital dis- 
tance would be approximately twice the distance at which 
the planet fills its Roche lobe (a ~ 2 ■ an, with au defined 
as the distance at which Rp = Rh, where Rh is the Hill 
radius). Using a Bayesian analysis of the sample of radial- 
velocity (i.e. non-transiting) planets, Ford & Rasio (20061 
found that the data was compatible with the closest plan- 
ets following a = a ■ an, the sini uncertainty on the mass 
of non-transiting planets preventing a determination of the 
numerical value of the multiplying factor a. 

On Fig. [3J the dotted line shows the locus of a = 2 • an- 
Clearly, this relation does not constitute a good descrip- 
tion of the present position of close-in planets on circular 
orbits. However, the slope of this relation seems to be a bet- 
ter description of the observed locus of circular orbits than 
the slope of the circularisation isochrone, with a ~ 3 — 4. 
One possibility to explain the larger proportionality factor 
in the relation between circularisation distance and an is 
that planets had much larger sizes earlier on. 

Intriguingly, the dependence on an seems to hold for 
hot-Neptune and super-Earth planets. These planets do not 
at all follow the mass-period relation of hot Jupiters, with 
the three presently known transiting super-Earth candidates 
all having extremely close orbits (P < 1 days), but once 
scaled to the planetary radius, their position align with a 
similar multiple of an ■ The data at this point are not suffi- 
cient for definite conclusion, and this can be tested as more 
small transiting planets are found and characterised. 



5 CONCLUSION AND PROSPECTS 

The empirical features above must be accounted for by 
any successful theory on the nature and history of hot 
Jupiters. The smooth relation between eccentricity and tidal 
timescale, as well as the pile-up of circular orbits at specific 
orbital distances at circularisation, strongly suggests that 
tidal circularisation and the stopping mechanism of close-in 
planets are related. The mass-period trend also suggests that 
heavier planets get circularised closer to the host stars, if at 
all. The lack of higher-mass planets on circularised orbits 
suggests that these planets end up in the star instead, pos- 
sibly because they provoke tides on the star that are strong 
enough to transfer angular momentum from the planet's or- 
bit to the stellar spin. 

The observed relation between orbital distance and in- 
trinsic parameters (size and mass) for planets on circular 
orbits shares some features with the predictions of |Ford ~&\ 
Rasio ( 2006 1 for inward migration caused by planet-planet 



scattering. 

The close association of migration and circularisa- 
tion suggest that the two mechanisms operate on similar 
timescales. If migration was much more rapid, tidal evolu- 
tion would not be efficient as a stopping mechanism. This 
argues in favour of secular dynamical interactions to explain 
the presence of close-in planets rather than migration in a 
protoplanetary disc. 

These are exciting times in the study of close-in exo- 
planets. As new objects are discovered by the dozen, the 
statistical inferences become more robust. Many features of 
the sample as a whole have become clearer: the abundance 
of misaligned systems and their relation with the nature 



of the primary, the role of inward migration, the relation 
between the energy budget and the planet size. These fea- 
tures seem to be related by an over-arching theme: most of 
the properties of close-in planets are actually controlled by 
their parent stars. When planets get closer than a few days 
in period, the star takes over their orbital properties, their 
atmospheric dynamics, and probably their size and lifetime 
as well. 
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